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LONG WRITE-UP 



1 Introduction 

Hypergeometric functions pFq appear in many branches of science. They ap- 
pear, in particular, in particle physics during the calculation of radiative cor- 
rections to scattering cross sections in loop [1-8] or phase space [9-11] in- 
tegrals. In the context of dimensional regularization, the parameters of the 
hypergeometric functions are usually functions of an arbitrary space-time di- 
mension D = 4 — 2e where e regulates infrared and/or ultraviolet divergences. 
For physical observables, only the limit e — > is of importance. Because di- 
vergences appear in form of poles l/e n , not only the constant term of the 
Taylor expansion around the physical dimension is of physical relevance. For 
this reason one is often confronted with the task of expanding hypergeometric 
functions about their parameters. 

Systematic approaches to the expansion of hypergeometric functions have been 
developed [12-14] and have been implemented in GiNaC [15,16], Mathemat- 
ica [17], and FORM [18]. These implementations, however, are restricted to 
the expansion of hypergeometric functions about integer-valued parameters. 

In computations involving massive particles [5, 19-30] the hypergeometric 
functions can contain half-integer parameters. Methods have been developed 
to expand hypergeometric functions with half-integer parameters [13,31-33]. 
However, none of these methods has been made available in a user-friendly 
package written in a widely used computer algebra programme. 

The aim of this work is to introduce a new method for expanding hypergeo- 
metric functions about half-integer parameters, and its implementation in the 
for this purpose extended version of the publicly available Mathematica [34] 
package HypExp [17]. 

This paper is articulated as follows. In section 2, which is quite technical, we 
introduce the algorithm that allows the expansion of hypergeometric functions 
about half-integer parameters. Section 3 describes the implementation of this 
algorithm in the existing Mathematica package HypExp. Section 4 contains 
examples of Feynman diagrams to which the package can be applied, and we 
conclude in section 5. Appendix A is devoted to the analytic continuation of 
the expansion, and appendix B contains the computation of the basis functions 
for several types of hypergeometric functions. 
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2 Expansion of hypergeometric functions about half-integer pa- 
rameters 

In this section we describe an algorithm to expand hypergeometric functions 
about half-integer parameters. The strategy of the algorithm will be to express 
the expansion of a hypergeometric function of a given type in terms of certain 
operators acting on the expansion of one hypergeometric function of the same 
type, the latter will be referred to as the basis function of this type. We 
will then deal with the task of expanding the basis function and carrying 
out the aforementioned operations. Our approach differs from other reduction 
algorithms [14,31-33,35-40] in that we relate only expansions to a given order 
with the expansion of the basis function, and not the functions themselves. 

We will start this section by giving some definitions and notation, before we 
explain the reduction to the basis function by means of an example, namely by 
considering a certain type of a 2 Pi function. We will then extend this method 
and show how to reduce a general function pF P _i to the respective basis 
function. Afterwards, we explain a strategy to find the all-order expansion of 
the basis functions and give some applications. 



2.1 Definitions and notation 

In this subsection we introduce some definitions and notation that will be used 
in the following sections. In this paper we consider hypergeometric functions 
(hereafter HF) 

pFp-i(Ai, A P ; Bi, Bp-i, z), 
whose parameters Aj and Bj have the following form 

1 

i + 7e or i + - + , 

where i is an integer, 7 is a real coefficient, and e is the parameter in which 
we will expand the HF. 

We will denote the finite part (i.e. either % or % + |) of a parameter by lower 
case Latin letters, the coefficients of the expansion parameter e will be labeled 
by Greek a or f3, and the parameter as a whole will be referred to by capital 
A or B. A, a, and a always correspond to an element of the first subset of 
parameters whereas B, b, and (3 correspond to the second subset of parameters. 

We will refer to an HF as being of the type Pj if i out of the P a fe 's and j out 
of the P — 1 frfc's are half-integers, and the other a^'s and bkS are integers. 
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We define the short-hand notation 

a:a a:b 6—1 



a:b 



a-1 



if a > b , 



j=a 



so that 



a:6 



r(6) = r(a) n(j) ■ 



(i) 

(2) 



jg.jg Algorithm 



The algorithm is based on the fact that one can write fractions of the type 



x' 



(i+j) n 

with i + j ^ in terms of the integration operator 





J + (j,n)[f](x) = (j + (j)) [J + (J>- 1)[/]] (*) 
in the following way 



(3) 



1 



dx'x'^x'^ J + (j)[y l ](x), 



— = J+(j,n)\y*\(x) 



(i+j) 1 

Polynomials in i and x can also be re-written using 

ix l = x-^-x\ i n x l = (x-^-\ x 1, , 
ax \ ax J 

for which we define the operators 



(4) 
(5) 



J !./) / !.r) ^^fir) EE J-(j,l)[f](x) , 
J-(j,n)\f\(x) = (j-(jj) [j-(j,n-l)[f]\ (x) . 



(6) 
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We then have 



ix* = J- (0) [y l ] (x) , i V = J" (0, ra) (x) . (7) 

We will use these facts to express the expansion of an HF in terms of these 
differentiation and integration operators acting on another HF of the same 
type, the latter will be referred to as the basis function of this type. In the 
following section, we will illustrate the algorithm with an HF of type 2\. 



2.2.1 Warm up: Hypergeometric functions of type 2\ 

We start from the definition of the hypergeometric function 2 F± 

*FMi,A 2 ,B u x)-l + r{AnA2) g r(Si + i)r(i + 1) * • (8) 

We consider the case where the finite parts of A 1 and B\ are half-integers and 
that of A 2 is an integer 

1 1 

A 1 — ai + - + a ± e, A 2 = a 2 + a 2 e, B 1 — b 1 + - + fi 1 e. 

Using the relation xF(x) — F(x + 1) we can transform the expression (8) into 



2-Fi(ai + ^ + aie, a 2 + a 2 e; b x + ^ + fte; x) 



r(| + ai£) r(i + a2£ )» ff( . + i + aie)lff( . + Q2£) 

3 3 

>:ai l:a,2 

j 1 j T(i + ± + aiejl (t + 1 + a 2 e) 



« t(* + i + HA0 r( 1 + | + /3lf) r (! + i) 1 



(9) 

where we made use of the short-hand notation defined in Eq. (1). The factor 
D in Eq. (9) can be turned into partial fractions with respect to i. This yields 
a sum of factors i n , + j + 7e) n or l/(i + j + 1/2 + ^e) n . 



6 



2 Fi(ai + ^ + «ie, a 2 + a 2 e; b x + ^ + /3ie; x) = 

=1+ n'o + |+ft £ ) r(4+Ae) 



n 1 / • ,i \ ^T~r 2 / ■ \ T(i + aie)r(l + a 2 e) 



oo 

X 



2: ( E , C k^n + E 



»=1 V>0,n v" 1 J 1 ' / j<0,n 

c 1/2 

r(» + l + ai ^ + | + a 2 e) , 

r(i + i + ^ lC )r(i + i) x ' 1 UJ 

where the coefficients Cf n and C j ra are polynomials in e. We will now show 
in turn for each of the four double sums how it can be expressed in terms of 
integration and differentiation operators acting on a basis function B. 

In the sum over positive j's the denominators can be expanded in e, which 
yields factors of + j) n that can be expressed as integration operators 

J+(j,n) = ^ fdx'x'i-^j (11) 

acting on 



B = r(| + /? ie ) ~ r(z + ^ + aie)r(i + l + q 2 e) f 

~r(i + a 1 e)r(l + a 2 e) ^ r(i + § + /3ie)r(i + 1) 

= 2 F 1 (^ + a 1 e,l + a 2 e,^ + (3 u x) -1 . (12) 

The same applies for the third sum since it contains terms of the form 
+ \ + j + 7e) n - They can be expressed using the integration operator 



x ' 



The factors i n in the last sum can be written as differentiation operators 

J-(0,n)=(x^J (13) 

acting on B. The second of the above sums requires more work. For these 
terms which have j = — k < we decompose the outer sum over i in two 
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parts, one running from 1 to k, and the other one running from k + 1 to 
infinity, according to the following structure, 



E 



X 



■ Bi = E p rz ^ + E • (14) 



^ (i - fc + 7 e) 
The last term can be worked out as follows 



■Bi 



E?-— w^+* = ^E 



7i (* + 7e 



i+k 



= ^ E W-w (-7^)' E 

1=1 V / 



=* fc E 

2=0 

oo 

=^ fc E 



Z!(n- 


1)! 


(n + l- 


-1)! 


l\(n- 


1)! 


(n + l- 


-1)! 


l\(n- 


1)! 



(- 7 e)' [j + (0,n + /)] ^x-V^j 
(- 7 e)' [j + (0,n + /)] (x~ k B(k,x)j , 



where 



0, i < k 
B^ i > k 



and 



(15) 



B(k,x) = '£x i B i , k . (16) 

i=i 

One can now write the terms with j < in D as an operator 



E r , w ^ - J+ U t)^ 



-E 



Ba l 



+ ^ E 



(n + Z-l)! 



^( i+i + 7e )n ^ Z!(n-1)! 



- 7 e) i J + (0,n + /)(^ J B(-j,x)) , 

(17) 



where -Bj is the i-th coefficient of the Taylor expansion in x of B about x = 
r(ai + \ + i + «ie)r(a 2 + 1 + i + a 2 e)r(&! + \ + fte) 



r(ai + | + aie)r(a 2 + 1 + a 2 e)r(oi + | + i + Pie)i\ 



(18) 



The final formula reads 
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2 Fi(ai + ^ + a i e , °2 + a 2 e, h + ^ + fte, x) 

0:6i 

n (J + \ + Ae) 

= 1 + 



X 



0:ai l:a2 

11 (j + \ + aie) II (j + « 2 e) 

Vi>°. n j<0,n,7 



+ E </ + (J + ^ n) + E C~ J- (0, n) )B . (19) 

It expresses the HF under consideration in terms of integration and differen- 
tiation operators acting on the basis function B which we have chosen to be 
as in Eq. (12). We will now generalize this procedure to arbitrary P F P _i. 

2.2.2 General method 

We now consider the general case for the expansion of 



P Fp_i(Ai, A P ; B u 5 P _i, x) = 

UL r(A t ) h T\U v{B 3 + i) r(< + 1) ' 1 Uj 



with 



Ai = a,i + - + oti€, 1 < i < r, Ai = di + c^e, r < i < P 

Bi = bt + - + foe, l<i<s, B i = b i + /3 i e, s<i<P-l, 



(21) 



where all Oj's and b^s are integers and the HF is therefore of type PJ according 
to section 2.1. Using the relations 

n:m 

xY{x) = F(x + 1), F(x + m) — T(x + n) J\(x + j) > 

i 

we can transform Eq. (20) into 



pFp_i(ai + i + «ie, ...,a P + a P e; b x + - + fte, ...,fc P _i + /5 P _ie,x) = 
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= 1 + 



( n r(m, + \ + Ae)) (^if r(m, + Ae)) 
(nr(m + | + aje)) ^_fi i r(n, + a,e)) 

n?=i n Oj + 1 + Ae) nf=-+i n U + Ae) 

x £i /_V £i / 

/ r nf.ai \ ( P nf.ai \ 

n n Ui + 1 + me) n n u + a, e ) 

V=l 3i ' / \«=r+l J; / 



oo 



I r ni'.ai \ / P n ( :a; 

n n (i + ji + 1 + «/e) n n (« + # + aze) 

Vi=l 31 ' J \l=r+l h 



X 



1 n?=i n (i+ji + i+Ae) nf= s +i n (*+^+Ae) 

D 
P 

. II 

(=1 " / Vi=r+1 J x l 



lir(i + ni + l + aie))l ft r(i + ra, + a,e) 

!=1 / \/=r+l 

(nr=i r(i + + 1 + Ae)) (n^i r (* + ™/ + Sj) + 1) ' 



(22) 

where the integers rij and m 8 can be chosen as appropriate since the choice of 
the basis function is not unique. The factor D can be partial-fractioned into 
a sum of factors i n , + j + 7e) n or + j + 1/2 + 7e) n . 

(23) 

Again, we write the terms in D as integration or differentiation operators 
acting on the basis function 



B r s ({n u ...,n P ;mi, ...,m P -i}, {a u ...,a P ;P u 
( jft r(m, + ± + Ae)) i T( mi + Ae)) 

fl !>, + I + ( n T(n ; + a,e) j 



1=1 / \l=r+l 



oo 



rilXi + raj + i + aje)) f _fi ^(i + n, + a,e) j ^ 



i (n?=i r(i + + § + Ae)) (nf= s +i r(i + + Ae)) r(i + 1) 
= pF P _x (Aq, iV P ; Mi, M P _i; x) - 1 (24) 

obtained by combining the first and last line of Eq. (22). To shorten the 
notation, we have used 
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Ni = rii + - + a^e, 1 < % <r, Ni = rii + ciie, r < i < P , 

Mi = wij + ^ + pie, 1 < i < s, Mi = mi + pie, s < i < P - 1 . 

(25) 

We see that the function B defined in (12) corresponds to B\({0, 1; 0}, {a±, cx<i\ Pi})- 
The final formula reads 



pFp-i (Ai, A P ; B>i,..., B P _i, x) 

n?=i n Ui + \ + Pie)) 1 1& n Ui + fc) 

31 / \ 31 y 

/ r nf.cn , \ / P nf.ai \ 

n n Ui + \ + a, e ) n non a, e ) 

v'=l Ji " / Vi=r+1 J! / 



= 1 + 



X 



+ E ^ + (j + i^) + E Cn (0, n) 

j,n Z n 

B r s ({ni, ...,n P ;m 1 , ...,m P _i}, ...,a P ;/3i, ...,/3 P _i}) 



(26) 



This formula relates a hypergeometric function to another one of the same 
type via differentiation and integration operators. It is therefore sufficient for 
the expansion of any pF P _i of type PJ to have 

a) the expansion of one HF of the type PJ (hereafter called the basis HF) 

b) a procedure to integrate and differentiate the expansion of this basis HF. 

The expansion of the basis HFs is treated in the next section. The various 
terms in the expansion will contain HPLs as well as rational functions. Since 
the HPLs are constructed by means of iterated integrations, they are well- 
suited for carrying out the required integration and differentiation procedures. 
We will therefore focus on the derivation of the expansion of the basis func- 
tions. 



2.3 All-order expansion of the basis functions 



The achievable depth of expansion of a hypergeometric function with half- 
integer parameters using our algorithm depends on the available depth of the 
expansion of the corresponding basis function B. In this section we describe 
one possible way to obtain an expression for the basis functions of some types 
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of HF to arbitrary order in the expansion parameter e. 



We first set up some notation. We introduce an associative, distributive but 
non- commutative product between HPLs of the same argument 

H(si, s n ;x) H(ti, ...,t m ;x) = H(si, s n ,ti, ...,t m ;x) 
(H(si... n ; x) + H(ti... m ] x)) H(ui...i] x) = H(si... n , x) + if (ii... m , u^.r, x) 
1 H(t u t m ; x) = H(t u t m ; x) , (27) 

with the short-hand notation Si... n = Si, s n . We can define a "division" with 
respect to this product 

H(s 1 , Sj, s n ; x) H(sj, s n ; x) = H(s u Sj_i; x) . (28) 

If the vector ti... m is not the last part of the vector s L .. n we define 

H( Sl ... n , x) if(fi... ro ; x) = . (29) 

The advantage of this notation is that the derivative "factorizes" 



^- (if (si... n ; x) if (ti... m ; x)) = ^if (si... n ; x)^J if(ti... m ; x) 
d j / \ 

— (if(si... n ; x) if(fi... m ; x)) = I — if(si... n ; x) J if (ti... m ; x) j < n . 

(30) 

For a j-fold derivative applied on a product whose left HPL has a weight 
smaller than j we have no factorization, but we have 

f^-(if(si... n ;x) 0if(t L .. m ;x))j 0H(t 1 ... m ;x) = l^-H( Sl ... n ;x) J . (31) 



£5.1 General strategy 
We make the ansatz 

(oo 
1 • » E c( Sl ,..., Sj ;x)// Sli ..., Sj .(/(x))] , (32) 
i=l si,...Sj=+,0,- 

with f(x) = \px for HFs of type P\ and 



v . / X 1 - y/l-X 

t(x)=tJ- or f(x) = ■= 

V; V 1-x l + v^ 



x 
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for HFs of type P£±i. These types of arguments have already been used pre- 
viously (see for example [13,23,26,31-33].) The function g(x) is given by the 
value of the HF with the expansion parameter e put to 0. We see that in the 
ansatz for HFs of the type P- ±1 we will get HPLs of imaginary arguments 
for real arguments of the HF to expand. The properties of HPLs of complex 
arguments are described in Refs. [26,41,42] as well as in appendix A. 

The problem reduces to finding the coefficients c aiv .. jS . ;a; in equation (32). Some 
of their properties can be stated. 

• The coefficients c slj ... jSj;x are homogeneous of order j in the a,, fa, that means 
that the powers of the different a's and fas sum up to the same power as e. 

• They have to be symmetric in a's and fas corresponding to equal a's or 6's. 

• Since an HF p + \F p reduces to an HF pFp_i if two of its parameters, one 
from the first and one from the second subset, are equal, we have conditions 
on the coefficients c Slj ... jS . ;x . If one of the a's, say Oj, is equal to one of the 
6's, say bj, the coefficient c should reduce to the coefficient of the reduced 
HF when we take the corresponding ctj and fa to be equal. 

The hypergeometric function pFp^i satisfies a differential equation T> of order 
P. Inserting the ansatz in the differential equation 

VB = 

the left-hand side can also be written as a sum of HPLs with coefficients 

oo 

E e 'E E V(s 1 ,...,s l )H Sl _ Sl (f(x)) = 0. (33) 

j=0 l s 1 ,...s l =+,0,- 

The differential equation is satisfied if all the coefficients V(s\, s n ) vanish. 
We consider the coefficient D of a given HPL with weight vector (vi, ...,v n ). 
Since the differential equation is of order P, this coefficient will get contri- 
butions from the coefficients c(...) of HPLs of weight n, n + 1, ...,n + P with 
weight vectors of the form (..., v i, ...,v n ). Using the notation (27), that is HPLs 
of the form, 

H( Sl , s k , f(x)) <g> H(y u s n , f(x)), < k < P. 

Thus, for the coefficient T>{y\, ...,v n ) we only need to consider a part of the 
ansatz B, namely 

B(vi, ...,v n ; x) = c(vi, v n )g(x) 

X ( 1 + E ei E ^s 1 ...f,x)H(s 1 ...ff(x))\®H(v l ,...,v n ,f(x)), 
V i =1 si,...sj J 

(34) 
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where we have defined 



c(s L . ,•; x) 



c(s L „ i ,77i, ...,v n ;x) 
c(v u ...,v n ;x) 



(35) 



and where the sum over the weights si, Sj runs over +, — , 0. The coefficient 
T>(vi, v n ; x) is now given by 



T>(vi, ...,v n ; x) = (VB(vi, v n ; xfj v n ; x) 



H(...)=0 



(36) 



The only dependence of the right hand side on the vector ...,f n ) is in the 
overall factor c(si, Sj; x) and, implicitly, in the coefficients c(si...j;x). Our 
strategy is now to construct coefficients c(si_f,x) such that the coefficients 
in Eq. (36) vanish. In the following, we illustrate this strategy by means of a 
simple example. 



2.3.2 Simple example 
We consider the HF 

2 Fi(i + «ie, 1 + a 2 e; ^ + /5ie;x) , (37) 

using the ± HPL weights defined in [41], we make the ansatz 

1 I 00 i 
B = -^~ U + E eJ E c(s l ,...,s j ;x)H[ Sl ,... lSj ;^)\. (38) 

We have taken g(x) = (1 — x)^ 1 since 

^(i.H !I ' = rb- (39) 

Since the expansion of Eq. (37) has to be real for negative values of x, we have 
to make sure that the coefficient of the HPLs in the expansion (which will 
have imaginary argument for x < 0) guarantees that the expansion is real for 
x < 0. We know (see [41]) that an HPL of imaginary argument is real if it has 
an even number of "+" weights and imaginary if the number of "+" weights 
is o Therefore we define 

y/xc Q (si, s n ) odd number of + in s l5 s n 
c e (si, s n ) even number of + in si, s n . 



1 provided the last element of the weight vector is not 0, which is not relevant here 
since such HPLs are divergent at x = and we are only interested in solutions of 
the differential equation regular at x = 0. 
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We have now for an even number of "+" in vi, v r , 

c e (v u ...,v n ) 



B(vi, ...,v n ;x) 



l-x 



x ^1 + e (yEc e (+)H{+; y/x) + c e {-)H{-- y/x) + c e (0)#(0; 

+ e 2 ^c e (+, +)#(+, +; y/x) + y/x~c e (+, -)#(+, -; y/x) 

+v^c e (+, 0)#(+, 0; V^) + y/x~c e (0, +)#(0, +; Vi) 
+c e (0, -)if (0, v^) + ce(0, 0)#(0, 0; Vx) 
+v^c e (-, +)#(-, +; v^) + Ce(-, -)#(-, -; v^) 

+c e (- 0)^(- ;v ^) M ®H{v 1 ,...,v n ; y /x) , (40) 

and 



5^ \ c (v u ...,v n ) 

B(v u ...,v n ; x) = 

1 — x 



X 



y[x~ + e (c (+)H(+; yfx~) + c (-)-v/a;#(-; y/x) + c o (0)y/xH(0; y/xj) 



+ e 2 ^c (+, +) v ^#(+, +; v^) + c D (+, -)#(+, -; v^) 



+c D (+, 0)if (+, 0; Vx) + c o (0, +)#(0, +; v^) + c o (0, (0, v^) 

+c o (0, 0)y/x~H(0, 0; v^) + c (-, +)#(-, +; v^) 

+ Co (-, -)yfcH(-, v^) + c c (-, 0)^(-, 0; v^)) j 

<g>if(ui,...,u n ; v^) , (41) 

for an odd number of "+" in vi, ...,v n . Inserting these expressions into Eq. (36) 
and equating the e coefficient of the left hand side to zero we get conditions 
on c e (+), c e (-), c e (0), c G (+) c (-) and c o (0), namely 



c e (+) = ai - /3i , c D (+) = a 2 

c e (-) = a 2 , c (-) = «i-/5i 
c e (0) = , c o (0) = -2ft . (42) 

These results are independent of the vector i> n , we can thus take this 

result to define a rule to construct the coefficients of the HPLs in our ansatz, 
Eq. (38) 



15 



c(+,W U ...,W n ) = 

a 2 odd number of + in {w±, ...,w n } 
cti — Pi even number of + in {wi, w n } 



c(wi, ...,w n ) x < 
c(0,w u ...,w r 

c(w u ...,w n ) 
c(-,wi, ...,w n ) 

c(wi, ...,w n ) x 



2 f3i odd number of + in {wi, ...,w n } 
even number of + in {wi, ...,w n } 

ai — Pi °dd number of + in {wi, w n } 
a 2 even number of + in {wi, ...,w n } 



(43) 



Using these rules we can show that 

V(v u ...,v n ;x) = 0, (44) 

for all vectors (vi,...,v n ). We can ensure that the boundary conditions are 
respected by setting 

c(0) = c(-) = ai c(+) = a 2 -Pi. (45) 

The first terms of the expansion in Eq. (37) are then 

1 1 

2 Fi(- + aie, 1 + a 2 e; - + Pie; x) = 

— ^— ( 1 + e (v^(«i - Pi)H(+; y/E) + a 2 H(-; y/x)\ 
+ e 2 ^a 2 («i - Pi)H(+, +; v^ 7 ) + >/x(ai - Pi)H(+, -; y/x) 



-2pi(ai - Pi)y/x~H(0, +; y/x) 



+ (ai - Pi) 2 ^H{-, +; y/x) + a 2 2 H(-, -; v^j j + C(e 3 



(46) 



This expansion could be rewritten in terms of the more usual HPLs of integer 
weights, yieling the same results as, for example in [13,32], but part of the 
structure would be lost. The success of the strategy for this case relies on the 
following facts 

a) The x dependence of the factors c(...;x) must be known and 

b) this dependence must be simple enough in order to minimize the vi, ...,v n 
dependence of the c(si. x) in Eq. (35). 



16 



2.4 Applications 



Following this strategy, we found explicit results for the all-order expansions 
of the following HFs. 



2 F 1 (^ + a 1 e,a 2 e;l + Pie;x) of type 2] , 

3 F 2 (i + a ± e, 1 + a 2 e, 1 + a 3 e; ^ + fae, 1 + (3 2 e; x) of type 3} , 

3 F 2 (a 1 e, 1 + a 2 e, 1 + a 3 e; ^ + foe, 1 + (3 2 e\ x) of type 3° , 



3 F 2 (- + aie, - + a 2 e, 1 + a 3 e; - + fte, - + /3 2 e; x) of type 3 



2 ' 



They are displayed in appendix B. Some of the expansions of those types listed 
at the beginning of the next section are too lengthy to be reported here. The 
result of the expansions found in all the expansions are expressed through 
HPLs of integer and ± weights of arguments 



y/x, 



Ik 



x 1 — yl — x 



1 — X 1 + y/l — X 



The analytic continuation of the expansion for these arguments is addressed 
in Appendix A. 

In section B.5 we show that the all order expansion of 

3-^2 Q + «ie, a 2 e, 1 + a 3 e; - + fre, ^ + f3 2 e; x^j of type 3 2 
cannot be expressed within the basis of HPLs used in this work. 



3 Extension of the Mathematica package HypExp 



We used the algorithm of the preceding section to extend the capabilities of the 
existing Mathematica package HypExp [17]. The second version of the package 
is now able to expand, in addition to HFs with only integer parameters, jFj_i- 
functions of the types 

r)2 r)l r)l r)0 

Z l> Z l> Z 0' Z l5 

33 o2 ql ol oO 

2 i °2> °15 °0' °1> 

4}, 4 

to arbitrary order in a small quantity about their parameters. The limitation to 
the number of classes of functions that can be expanded is due to the fact that 
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the all-order expansion of a basis function for each class has to be computed 
on a cases-by-case basis. The package is publicly available from [43]. The re- 
sults are displayed in terms of rational functions, logarithms, polylogarithms, 
Nielsen polylogarithms, and HPLs. To treat the harmonic polylogarithms, the 
package HypExp uses the package HPL [44] whose second version [41] imple- 
ments the harmonic polylogarithms of complex arguments, as required by the 
expansion about half-integer parameters. The results given by the package 
are valid in the entire complex plane. However, along the branch cut running 
along the real axis from 1 to +oo care has to be taken in order to reproduce 
the correct imaginary part, see section 4 as well as appendix A. 

In this section we shortly describe the new features of the package HypExp, for 
a more detailed description, we refer to our previous work [17]. 

After installation^, the package HypExp may be loaded using the following 
command, 



<< HypExp'' 



This should be done at the beginning of the session. 
HypExp [Hypergeometric2Fl [ . . . ,x] ,e,n] or 

HypExp [HypergeometricPFQ [ . . . ,x] ,e,n] returns the e-expansion of the en- 
closed hypergeometric function (HF) through order 0(e n ) if the type of HF is 
supported. 



HypExp [Hypergeometric2Fl [— + e, 1, 2-e, x], e, l] 



■ 2 e 



/ , I r 1-Vl-X i r 1 - V 1 - X 

3 -1+Vl-x + x -4 (-1 + x) Log 1 +2 (-1 + x) Log 1 + 

1+Vl-x 1 + Vl-x 



HypExp [ Hypergeometr icPFQ [ j— + e, -l + 3e, lj, j e, 1 + eJ, x] , e, l] 



1 + e 4 x (11 - 9 Log [1 - x] ) - 3 Log[l - x] 

3 .fV 9 



The function HypExp applied to anything else but a HF will return the same 
object. The result is not returned as a SeriesData object because this would 
have the effect of forcing the expansion of the rest of the expression. 

The prefactors that accompany the variable e can also be symbolic, and the 
expansion also works for argument x = 1 as shown by the following examples, 



see link [43] for more information on the installation procedure 
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HypExp [HypergeametricPFQ^j— + a e, ltjSe, lj, j— +re, l + <5ej, xj , e, lj // FullSimplify 

— - — -2yeHPL[{0, plus}, Vx~l + 
2 V^T 1 



(a - /3 - y + 5) e HPL {minus, plus} , V~x~ + Log -1 - 



-1 + Vx~ 



(1 + 2 ye + (-13 + 6) e Log[l - x] ) 



HypExp [HypergeametricPFQ^jl + e, 2 e, 3eJ, j— +2e, ej, lj , e, 2] 



1 


19 Log[2] 


) + ^ 


'161 


5/T 2 


19 Log[2] 


Log [2] 2 ' 


— + e 












6 


60 3 






9 


30 


3 t 



The package also updates Series to allow it to expand compound expressions 
containing hypergeometric functions. The difference between this and the func- 
tion HypExp is that the other functions of e are also expanded. 



Series ^Hypergeometric2Fl J — + e , , e, xj, {e, 0, 1}] 



i Log 



2 Vx 



4 Vx~ 



+ 2 x 

V -1 + x V -1 + X 



■ 2 i HPL {0, plus}, i . 



1 - x 



- 3 i HPL {minus, plus } , i 



1 - x 



-1 + x 



Log 1 + 



■ 4 x 



1 - x J 



1 + i • 



Log 1 + + i Log 

-1 + x L 1 - x 1 1 



1 - i ■ 



1 + i 



1 + i ■ 



4 i x Log 



- 2 i Log 1 + 



1 - i 



1 - x J 



Log 



1 - i 



e + 0[e] 



Series [— HypergeometricPFQ [ j — + e, 1, l-2ej, jl - 3 e, ej, xj , {e, 0, 0}j 



ArcTanh Vx 



2 Vx~ 



:hpl[{0, plus}, Vx~] + HPL [ {minus, plus}, Vx~] - 2 Log [ - — -^^j -Log[^ — 



Log [1 - x] 



1 - Vx 



" 1 - Vx 



■Ole] 1 



The use of the function HypExp is preferable in particular if the expansion or- 
der is high or the expression containing the HFs is large. The results given 
by the package are not systematically simplified using Simplify, since the 
simplification might take longer than the expansion itself, in particular for 
expansions to high orders. A Simplify might produce a more compact result. 
For additional commands and functions implemented in the package see [17]. 
The package was developed in Mathematica 5 and tested on both versions 5 
and 6. 
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S = (Pi + p 2 



2 2 

pj = m 




2 2 
£>2 = W 



pl = o 



4~ 



q 2 = (pi + pa) 



p| = 



Fig. 1. Left panel: One- loop vertex correction with one massless line (dashed) and 
two massive ones (solid) with equal masses. Right panel: Three-loop master integral 
with six internal massless lines. 

4 Examples and applications 



A first application of the present package has already been given in Ref. [29] 
for the case of a three-loop master integral with HQET propagators. Four 
more examples with up to four loops are given below. 



4-1 One-loop vertex correction 



Our first example is the one-loop vertex correction displayed on the left in 
Fig. 1, of which a one-dimensional Mellin-Barnes representation was derived 
in Ref. [45]. Assuming all propagators to be raised to unit power, we get 



[dk] 



1 



[k 2 — m 2 ] [{k — Pi — P2) 2 — rn 2 \ [(k — Pi) 



2e 



dz T(-z ) T 2 
2tH 



-e-z)T(l + e + z) 



T(-2e-2z) 



m 



[s - irj[- l - e - z (47) 



with 



[dk] 



D/2 



d°k 



(4£ 



and D = 4 — 2e. Here and in the following we tacitly assume that all propaga- 
tors contain an infinitesimal +irj. The integration contour in Eq. (47) is chosen 
such that it separates left poles of the T-functions from right ones. Applying 
the formula 



T(2z) = ^T(z)T(z+ 1 - 

WIT I 



(49) 
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to the T-function in the denominator of Eq. (47) and subsequently summing 
the residues to the left of the contour, we obtain 



/ = ir r ( i + e) 



m 



-1-6 , 3 x 

2^(1,1 + 6; -; s), 



(50) 



with s = (s + irj) / (4m 2 ) . We demonstrate below how the e-expansion of this 
expression can be performed with the package (we set m— 1). 



vertexcorrection = Exp [e EulerGamma] Gamma[l + e] / 2 / e Hypergeometric2Fl [1, 1 + 6, 3/2, s] 

e EulerGammaiS Gamma [e + 1] Hypergeometric2Fl [ 1, e + 1, -§-, s] 
2e 



Timing [HypExp[#, 6, 3] & //©vertexcorrection] 

{2 . 6216 Second, 

e EulerGamma 6 Gamma [ £ + X ] 



1 

2 6 



1 



1 - s 



HPL[{minus, plus}, i / - S - ] 



Log [ — — — + l] + 2 HPL [ {minus, plus, minus} , i / — — — ] Log [ — — — + l] 



2 HPL [ {minus, plus, minus, minus} , i / — — — ] e 



-HPL [ {minus, plus , minus } , i / — - — ] - 



HPL [ {minus, plus } , i / j S ^ _ ] Log [ - S - + l] 



HPL [{minus, plus}, i / - S - ] 6 



Log 



Series [%[ [2]], {s, 0, 2}] //Simplify 

L °g[- j.y^j^] i^-^r HPL[{minus, plus}, i ^ - ] 



1 

48 s 



4 s E 



-1 + s 



4 s 



6HPL[{minus, plus, minus}, i / - ] 



6HPL[{minus, plus}, i I - ] Logfy^] - n 2 Log[- f =H 



1 

48 s 



: HPL I {minus, plus, minus, minus}, i 



4 HPL [ {minus, plus, minus}, i / ] Log [ — - — ] + HPL [ { minus , plus}, 



-1 + s 



|/T 2 + 2 Log[ - \ ] | - 2 Log[ y ~_ 1+s j - ..[:., 



+ 0[e] J 
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The result of the expansion agrees with the one recently given in Ref . [46] . By 
means of the commands HPLpm21ml and HPLConvertToKnownFunctions from 
the HPL package [41,44], the above expansion can be converted to more com- 
mon functions like logarithms, polylogarithms, and Nielsen polylogarithms. As 
a cautionary remark, we would like to stress that above threshold the sign of 
the ^-prescription in Eq. (50) is crucial in order to obtain the correct result, 
see appendix A and Refs. [41,44]. 

4-2 Massless three-loop master integral 

Our second example is the computation of the massless three-loop master 
integral which was identified in Ref. [8] and which is depicted on the right in 
Fig. 1. After Feynman parametrization, this integral reads 




1 



(51) 



1 








x ZT 3e (l - r)~ 3e (l - t)~ e (l - y) e {l - zf- 2e r- e z 
x ll+r((l- z)zD 2 + (1 - y )(l - z )tD 




-3e 



(52) 



with 



D = 1 




(53) 



After numerous variable changes we get 




A/~2 8£ - 2 7rr 2 (l - 3e)r 5 (l - e)r(3e) 
e r(2-4 e )P(§-2e) 




(54) 
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We see that we have an HF of type 3° with one half integer parameter. We 
can expand each factor in e. The powers of s and (1 — s) can be expressed 
using HPLs (and HPL weight functions) 



{1-8) 



)- 3e = ]T3ViJ(n;s) 



3 



= 1 + 3eif(l; s) + 9e 2 H(l, 1; s) + 27e 3 #(l, 1, 1; s) + . . . 
(sr = £(-l)V#Cl; S ) 



= 1 + 2eH(0; s) + Ae 2 H(0, 0; s) + 8e 3 # (0, 0, 0; s) + . . . 



-l+2e _ 




(55) 



The HFs can be expanded with HypExp. The first one includes only integer 
parameters and hence its expansion will contain only HPLs of argument s. 
The second one has half-integer parameters and will have HPLs of square root 
arguments 



so that we can convert the HPLs of this argument into HPLs of argument s 
by applying (twice) the command HPLConvertToSimplerArguments from the 
HPL package. The next step is to expand the product of HPLs into a sum of 
HPLs which we can then integrate by means of the integration routines of 
HPL. This procedure is not restricted to a specific depth of the expansion, so 
we could, in principle, expand A & 2 to all orders. After factoring out an ap- 
propriate combination of prefactors, the result can be written as the following 
homogeneous sum, where the power of e grows with the transcendentality of 
the coefficient. 




with 



x = 




Since s is in the interval (0, 1), we have 
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u,z (l-5e)(l-4e)e 
7tt 



x 



2C(3)- e ^ + e 2 ( VC(3) - 10C(5)) 

+e 5 (-^C(5, 3) + 1000C(3)C(5) - 22vr 2 C(3) 2 + + 0{^) 

(56) 

where we have encountered a multiple zeta value in the last term. 



4-3 Two-loop vacuum diagram with three massive lines of equal masses 



Our third expample is the two-loop vacuum diagram with three massive lines 
of equal masses which is displayed on the left in Fig. 2. It has been considered 
at several places in the literature, see e.g. Refs. [19,25,27,47-50]. In Ref. [19] 
the following expression for positive integer values of ni, n 2 , and n 3 was derived 
in terms of hypergeometric functions, 



al D k ld u l 



1D1 



[k 2 - m 2 ]" 1 [P - m 2 ]™ 2 [(k - I) 2 - m 2 ]" 3 



,l-ni23 ( m ^D-n 12 3 



f T(D/2 - n 3 ) T(n 13 - D/2) T(n 23 - D/2) r(n 123 - D) 

x I t(d/2) r(m) r(n 2 ) i> 1233 - d) 

n 3 , n 13 - D/2, n 23 - D/2, n l23 - D 
(n 1233 - D)/2, (n 1233 - D + l)/2, n 3 -D/2 + 1 
T(n 3 - D/2) T(n 12 - D/2) 



x 4 F 3 



r(n 3 ) T(n u ) 



( 



x 4 F 3 



n u n 2 ,n 12 -D/2,D/2 



\ (m 2 )/2, (m 2 + l)/2, D/2 - n 3 + 1 
with the usual abbreviation 



(57) 



ni + ... + n,- 



(58) 



It turns out that for positive integer values of rii, n 2 , and n 3 the 4 F 3 functions 
in Eq. (57) can be expressed as a linear combination of 2 F 1 functions of type 
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2^, which can then be expanded by means of the package. For instance, for 
n i = n 2 = n 3 = 1 we obtain the known result [19,25,27,47-50] 



2e 2 

+4*(-L3,!)-^-^(3)) + 0(aj, m 



with 



In Eq. (59) we have used a slightly different notation compared to the results 
given previously in the literature [19,27,47-50]. Relations among the different 
notations involve (see also Refs. [25,51]) 



^(5) = ^ + ! S 2 = ^ + 3V3Ls 2 (f) 

= |! + ^l.n(3) + f»(-iM) <«) 

! jf* cfo m|2sin(§)| (63) 



So 



9^3 

Lsj(0) = -jf dx ln^ 1 |2sin(f)| (64) 



$(Z,8,a) = £ [(fc + a) 2 K 2 - (65) 



4-4 Four-loop tadpole with three massive lines 



Our last example is related to the previous one, namely the four-loop tadpole 
with three equal massive lines displayed on the right in Fig. 2. We follow again 
Ref. [45], where the following one-dimensional Mellin-Barnes representation 
was derived 
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Fig. 2. Left panel: Two-loop vacuum diagram with three massive lines of equal 
masses. Right panel: Four-loop tadpole diagram with two massless lines (dashed) 
and three massive ones (solid) with equal masses. 



T{n u n 2 , n 3 , n 4 , n 5 ) = J [dfei] J [dk 2 ] J [dk 3 ] J [dk A ] 

1 

X Wi + k 2 ) 2 - m 2 T 2 Kh + h + h) 2 - [kt} n4 [k 2 3 ] n5 

(_l) n l2345 ^ m 2^8-4e-ni2345 g 4e7B 

= r( ni )r(n 2 )r(n 3 )r( n4 )r(n 5 ) 

4-i oo 

r(2-e-» 4 )r(2- e -n 5 ) ,iz_ , 
r(2-e) J 2ni v ; 

— i oo 

x r(-«)r(2 - e - n 2 - z)T(A - 2e - m 2 - 2) 

r(-6 + 3e + n 1245 + z)r(-8 + 4e + n 12345 + 2) .... 
X T(4-2e-n 12 -2z) ' (66) 

with the same integration measure as in Eq. (48). We shall again restrict 
ourselves to the case in which all the rij as well as the mass are equal to unity. 
We apply again formula (49) and obtain 



e 4e ^r 2 (l-e) 2f . 
7X1,1,1,1,1) = ^ J 2 2 - « K 

+ i 00 

' dz 



j^-T{-z)T{2-z-2e) 

— i 00 

F{l-z- e)T(z + 3e - 2)T(z + 4e - 3) (I 



r(|-z-c) 



(67) 



This expression can now be converted to a Meijer-G function [52-56] and sub- 
sequently be expanded in terms of hypergeometric functions. This procedure 
is in this case equivalent to summing all residues of left poles of T-functions 
in Eq. (67). The result can be displayed in the following closed form 
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1X1,1,1,1,1): 



2 3-4 £ e ^ E7rT 2 (l-e) 
sin(7re)r(2 - e) 



V^r(e)r(-l + 2e)r(-2 + 3e) 



r(2- e )r( 



26) 



x 3 F 2 (e, -1 + 2e, -2 + 3e ; 2 - e, -\ + 2e ; J) 
r(-| + e)r(-2 + 3e)r(-3 + 4e) 
r(-f + 3e) 

x 3 F 2 (-l + 2e, -2 + 3e, -3 + 4e ; e, -| + 3e ; £ 



(68) 



which is well-suited for expansion with the HypExp package. After some ma- 
nipulations and simplifications on the harmonic polylogarithms, one gets for 
the expansion in e up to the finite part 



1 1 / 97 7T 

T(l, 1, 1, 1, 1) = -f- + ^ + — + 

v ' ' ' ' > 4e 4 e 3 \48 



2\ ^ 

12) ? 
1 



977T 2 4C 3 7T 4 
+ — + — + 

144 3 12 1728 



833 
288 



7T 2 ( 3 \ 1 4177 
3 3 / e 432 



99 + 16tt z 



24^W(|) 



+ 0{e) , 
(69) 



in agreement with the findings of Refs. [45, 57, 58jj0- As a closing remark, we 
mention that hypergeometric functions of argument z = 1/4 as in our last two 
examples have been studied extensively in the literature and the structure of 
the coefficients in their expansion has been analysed [25,59,60]. 



5 Conclusion and outlook 



In this paper we have presented a new algorithm to expand hypergeometric 
functions about half-integer parameters. The strategy is to express a HF of 
a given type in terms of integration and differentiation operators acting on a 
particular HF of the same type, the latter is called the basis function. The 
choice of the basis function for a given type of HF is not unique and can be 
taken such as to be as convenient as possible. The method therefore requires 
as ingredients the knowledge of 

a) The expansion of one HF of the same type as the one to be expanded, 
namely the basis function. 

b) Integration and differentiation routines for the functions occurring in this 
expansion. 

The method is quite general and might be applied to HFs p-Fp_i of higher P 
or to other types of HFs. 

3 Note that s 2 in Ref. [58] is defined differently from 5 2 in Eq. (63) 
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We explain a strategy to compute explicit all-order expansions of the hyper- 
geometric functions and give explicit results for 

+ <xie, a 2 e; 1 + fiie;x) of type 2j , 
3 F 2 (i + a ± e, 1 + a 2 e, 1 + a 3 e; ^ + P ± e, 1 + /3 2 e; x) of type 3} , 

3 F 2 (o;ie, 1 + a 2 e, 1 + «3e; ~ + Pic, 1 + /^e; x) of type 3? , 

11 11 

3 F 2 (- + one, - + a 2 e, 1 + a 3 e; - + fte, - + (3 2 e; x) of type 3 2 . 

By taking appropriate limits, these expansions yield the expansion of the HFs 



1 1 
2 Fi(- + «ie, 1 + a 2 e; - + fre; x) 

2 Fi(l + ^e, 1 + a 2 e; ^ + fte; x) . 

The algorithm described here has been implemented for HFs of the types 

r\2 r)l r\l r)0 

Z l5 Z l5 Z 0> Z l' 

q3 q2 ol ol oO 

°25 °2' °1? °0> °1> 

41, 4^ 

in the already existing Mathematica Package HypExp. 

In order to show the relevance of hypergeometric functions with half-integer 
parameters and to demonstrate the performance of the package we have given 
examples of Feynman diagrams with up to four loops. We reproduced known 
results, but also gave new results as for instance in the case of the massless 
three- loop master integral A 6t2 . 
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A Analytic continuation of the expansion 



The results given by the package are for arguments in the interval (0, 1). When 
the expansion has to be evaluated elsewhere in the complex plane, care has 
to be taken. Let us analyze the three possible types of arguments inside the 
HPLs. 



A.l H(...,iJ^r a 



The expansion of HFs of the types 2°, 3°, 2f, and 3^ contains terms of the 
form 



H ...,z< 



x 



1 — X 



as, for example in 



.2 '2' 2' 
The function 





in the argument of the HPL maps the complex plane into the upper half plane, 
with a branch cut starting at z = 1 and following the real axis to +oo, as does 
the hypergeometric function. For x < 0, the argument of the HPL is real, but 
there is a sign ambiguity due to the convention one uses to define the square 
root of a negative number. Since HPLs with an odd number of + weights are 
odd, this ambiguity is removed by the factor of y/x accompanying the HPL, 
provided the same convention for the sign of the square root is used. 



A. 2 H(...,y/x) 



The expansion of HFs of the types 2\, 3\, 3|, 4\, and 4| contains terms of the 
form 

#(..., V*)- 

for x < 0, the argument of the HPL is complex. HPLs of complex arguments 
with ± weights are either purely imaginary or purely rea{3 depending on the 
number of + weights in the index vector. The HPLs that are purely complex 
are accompanied by factors of \fx and the HPLs that are real are not, so that 

4 provided there is no weight on the right of the index vector, which is the case 
in the expansion we deal with. 
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the expansion is also real for x < and, again, the sign ambiguity is removed 
if the same convention for the square root of a negative number is used both 
in the argument of the HPL and in its prefactor. 



The expansion of HFs of the types 2j and 3j contains terms of the form 



In this case, the HPLs are more conveniently expressed with integer weights. 
The branch cut of the argument of the HPLs for x > 1 corresponds to that of 
the hypergeometric function. There the argument of the HPL is on the unit 
circle and the HPL develops an imaginary part, as does the HF. 



B Computation of the basis functions 



Since the integration and differentiation operators introduced above with ap- 
propriate coefficients allow one to go from one HF to another of the same type 
Pg, the choice of which HF should be the basis function is arbitrary. In this 
section, we show the result of the application of the strategy described in sec- 
tion 2.3.1 to functions pF P _ 1 of types 2j, 3|, 3?, and 3%. We will finally show 
in section B.5 that in the expansion of HFs of type 3^ new weights have to be 
introduced. Some of the expansions in this section could have been also com- 
puted with other methods like for example those presented in Refs. [13,31,32], 
but our method yields the expansion directly in a form suitable to be used in 
the algorithm described above. 



B.l HFs of type 2* 

For the expansion of the basis function £>o({0, 0; 1}, {a±, a 2 ] Pi}) we consider 
the expansion of the HF 

2F1 Q + «ie, a 2 e; 1 + fre; x^j , (B.l) 

for which we found the expansion 

/l \ 00 

2 Fi (- + a 1 e,a 2 e;l + p 1 e;x =l + £^ ]T c(l)H t (u) . (B.2) 

V y j'=i ze{-i,o,+ip' 
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The coefficient functions c are found recursively by applying the following 
rules. 



c(0) 


= c(l) = 0, 


(B.3) 


c(-l) 


= 2a 2 , 


(B.4) 


c(l,...) 


= 0, 


(B.5) 


c(-l,...) 


= 2a 2 c(...) , 


(B.6) 


c(0,-l,...) 


= 2a 2 (ai - a 2 )c(...) + 2aic(0, ...) , 


(B.7) 


c(0,l,...) 


= 2(ai + a 2 - fo) (a 2 c(...) + c(0, ...)) , 


(B.8) 


c(0,0,...) 


= -Ac(o,...), 


(B.9) 



where 



LU 



1-vT 



1 + a/T" 



(B.10) 



5.5 HFs of typeS] 



We consider the function 



3-^2(7; + "if, 1 + « 2 e, 1 + a 3 e; ^ + fte, 1 + /3 2 e; x) . 
The expansion of this function is best written as a function of the parameters 



s = a 2 + a 3 — fo , 
d\ = a 2 — (3 2 , 
d 2 = a 3 - fo , 
d 3 = ai - fo. 

Like in the example in section 2.3.2 we make the ansatz 



(B.ll) 



3-^2(7; + «ie, 1 + a 2 e, 1 + a 3 e; ^ + foe, 1 + /3 2 e; x) 



l + E eJ E c(Z;x)!T,(V5) 



1 -x 



(B.12) 



j'=i /e{+,o,-}J 

The coefficient functions are found recursively by applying the following rules 



c(l; x) —c(l) x < 



ix, odd number of + in I 
1 even number of + in / 

c(+) = d 3 , c(0) = 0, c(-) = s. 



(B.13) 
(B.14) 
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The prefactor for + weights makes sure that the expansion is real, as adding 
a "+" in the weight changes a real HPL into an imaginary one. 



c(+,w u ...,w n ) = 

s odd number of + in {w±, w n } 



c(wi, ...,w n ) 
c(0,wi, ...,w r . 



c(w u ...,w n ) < 



C{-,Wi, ...,w r , 



c(w u ...,w n ) < 



d 3 even number of + in {w 1 , ...,w n } 



odd number of + in {w±, ...,w n } 



—2ft even number of + in {wi, ...,w n } and wi = 
— — - even number of + in {wi, ...,w n } and wi ^ 



d 3 odd number of + in {wi, ...,w n } 
s even number of + in {wi, ...,w n } . 



(B.15) 



The basis function is then found by subtracting unity from the HF. 



^({0, 1, 1; 0, 1}, {ai, a 2 , a 3 ; ft, ft}) = — - [x + E ^ E ^)^(v^) 

X V i =1 ie{+,o,-}' 

(B.16) 

The expansion of the basis function B\({0, 1; 0}, {cti, a 2 ; ft}) is recovered by 
setting 

s — > a 2 , ^2 — > 0, rf 3 — > «i - ft, 
and that of the basis function -B°({1, 1; 1}, {cti, a 2 ; ft}) by taking the limit 



s — > cti + a 2 — ft, di — > ai — ft, <i 2 — * a 2 — ft, d 3 







We see that in this limit all HPLs with weight + disappear from the expan- 
sion. Since only HPLs with weight and — can be transformed into HPLs of 
argument x 2 (see [41]) we expect only such HPLs to be present in the expan- 
sion of HFs with only integer parameters, since we know that the expansion 
of these HFs can be written in terms of HPLs of argument x [14]. 



B.3 HFs of the type 3? 



We consider the HF 3 F 2 (otie, 1 + a 2 e, 1 + a 3 e; | + fte, 1 + fte; x) and make the 
ansatz 



32 



3 F 2 (a 1 e, 1 + a 2 e, 1 + a 3 e; - + fte, 1 + (3 2 e; x) = 



1 + E <l;x)H t 



x 



1 — x 



(B.17) 



3=1 ie{+,Q,-Y 

The coefficient functions are found recursively by applying the following rules 



i^Jjz^, °dd number of + in / 
1 even number of + in / 
c(+) = -a 3 , c(0) = 0, c(-) = 0, 



c(Z; x) = c(l) x < 

C ( + ) = -a 3 , C ( 

c(+,wi, ...,«;„) = c(wi, ...,«;„) 



(B.18) 
(B.19) 



x < 



d x d 2 - s(d 3 + (3i) 



odd number of + in {wi, w n } 



even number of + in {wi, ...,w n } and W\ = + 



s — d± — d 2 — ds + j3\ even number of + in {w\, w n } and w\ ^ + 
c(0,wi, ...,w n ) = c(w 1 , ...,w n ) 

odd number of + in {w±, ...,w n } 
even number of + in {w\, ...,w n } and w\ = + 



x < 



-2fc 
2d x d 2 



2(di + d 2 — s) even number of + in {w±, w n } and W\ ^ + 
c(-,w x , ...,tw n ) = c(wi, ...,«;„) 

-s — a 3 odd number of + in {w 1 , ...,w n } 

< h-l even number of + in {w±, w n } and W\ = + 



x < 



d\ + d 2 — s even number of + in {w±, w n } and w± ^ + , 



where we used the definitions of (B.ll). The basis function 
B° 1 ({0,l,l;0,l},{a 1 ,a 2 ,a 3 ;f3 1 ,p 2 }) 

is then simply 



(B.20) 



5°({0, 1, 1; 0, 1}, {ax, a 2 , a 3 ; fa, &}) = ]T e J £ c(Z; x)^ 



j=i ie{+,o,-}' 



x — 1 



(B.21) 

The expression for £?°({0, 1; 0}, {«!, a 2 ; is found by setting a 3 — ► /3 2 in 
(B.21). 



33 



B.4 HFs of type 3j 



We consider the HF 3 F 2 (| + «ie, \ + a 2 e, 1 + a 3 e; \ + (3±e, \ + /3 2 e; x) and make 
the ansatz 



11 11. 

3 F 2 (- + «ie, - + a 2 e, 1 + a 3 e; - + foe, - + /3 2 e; x) 

oo 



1 x V i =1 ie{+,o,-}i V V 1 



x 



X 



(B.22) 



The expansion of this function is best written as a function of the parameters 
s — cci + a 2 - (3i - fa , 

S/3 = Pi + @2 , 

p a = a>ia 2 , 

Pp = Pifa. (B.23) 
The coefficient functions are found recursively by applying the following rules 



c(l;x) =c(l) x < 



*yT^r ' number of + in / 

1 even number of + in / 
c(+) = s, c(0)=0, c(-)=0 



(B.24) 



c(+,wi, ...,w n ) = 

!a 3 odd number of + in {wi, ...,w n } 
s even number of + in {w±, w n } 

c(0,wi, ...,w n ) = 

—2spc(wi, ...,w n ) — 4p ( gc(w 2 , w n ) odd number of + in {wi, w n } 

< 

even number of + in \wi, ...,w n } 

c(-,wi, ...,w n ) = c(w x , ...,w n ) 

s odd number of + in {wi, ...,w n } 
< ' (B.25 

a 3 even number of + in \wi, ...,w n } . 
The basis function £?f ({0, 0, 1; 0, 0}, {cci, a 2 , a 3 ; (3i,{3 2 )} is then simply 
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^({0,0,l;0,0},{ai,a 2 ,a 3 ;/3i,/fe}) 



The results of Eq. (43) are recovered by setting either cti or a 2 equal to one 
of the fas. 



5.5 HFs of type 3\ 



For the expansion of the basis function ^({0, 0, 1; 0, 0}, {cci, ck 2 , 0:3; fa}) 
we consider the expansion of the HF 

1 11 

3 F 2 (- + «ie, a 2 e, 1 + ot 3 e\ - + fae, - + fae\ x) , (B.27) 

for which we found the expansion 



(\ 1 1 

zF<2 V2 + " ie ' a2C ' 1 + a36 ' 2 + 2 + ^ 26 ' X 
1 — a^tioH (+; a;) 

- «2e 2 ((- («2 + "3) vH(—, +; a;) + a 3 H(+, +; cj) 
+2 - ft -/^wi^O, +;«;)) 

+ a 2 e 3 ^a 2 a 3 ujH (+, +, +; a;) - (a 2 + «3) 2 uH(—, -, +; a;) 

-4 (ai - /3i - /3 2 ) 2 wif (0, 0, +; u) + a 3 (a 2 + a 3 ) #(+, -, +! 
+2 (a 2 + a 3 ) («x - ft - ft) wif (0, +; u) 
-2a 3 (ai - A - ft) #(+, 0, +; u) 

+2 (a 2 + a 3 ) ( ai - ft - ft) wif (-, 0, +; w)^ 
+ 4e 3 a 2 ( ai - A) (ax - ft) uoH Q, 1 +; + 0(e 4 ) , (B.28) 
where 



UJ = Ik 



, 1-x 

In the last term, we have to define new weights for the HPLs: 

/i(t) = — f^=. (B.29) 

The appearance of this new weight shows that the HPLs of weights +, — , are 
not sufficient for the construction of the expansion of all HFs with half-integer 
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parameters. This is not a limitation to the algorithm described in this paper. 
However it is a limitation for the current implemetation, as these weights are 
not supported by the package HPL. 
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